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We study abrupt changes in the dynamics and/or steady state of fermionic dissipative systems 
produced by small changes of the system parameters. Specifically, we consider open fermionic 
systems whose dynamics is described by master equations that are quadratic (and, under certain 
conditions, quartic) in creation and annihilation operators. We analyze both phase transitions in 
steady state, as well as "dynamical transitions" . The latter are characterized by abrupt changes in 
the rate at which the system asymptotically approaches the steady state. We illustrate our general 
findings with relevant examples of fermionic (and, equivalently, spin) systems, and show that they 
can be realized in ion chains. 



I. INTRODUCTION 

Motivated by the impressive experimental control over 
many-body quantum states and dynamics open 
many-body quantum systems have received increasing 
experimental and theoretical attention in recent years. 
On the one hand, the decoherence introduced by cou- 
pling to an environment is a major challenge to quantum 
information processing Q, on the other hand, it can play 
a constructive role for quantum computing a, |4| , state 
preparation [5H8I1 , entanglement generation J8 ■jlCM . quan- 
tum memories or quantum simulation [12Hl6j . 

These exciting possibilities drive the interest in under- 
standing the steady state phase diagram of open systems 
in detail (l7| . Of particular interest are points of transi- 
tions between different phases of the system. For closed 
systems at zero temperature, the phase diagram and 
quantum phase transition can be understood by studying 
the low-lying energy eigenstates of the system's Hamil- 
tonian [18]. In particular, the non-analyticity of certain 
expectation values as a function of an external parame- 
ter, that characterizes the quantum phase transition, can 
only occur if the gap of the Hamiltonian closes, i.e., the 
energy difference between ground state and first excited 
state vanishes. Quantum phase transitions are thus de- 
termined by the low energy spectrum of the Hamiltonian 
governing the dynamics of wave functions 
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In this paper we study abrupt changes in the phys- 
ical properties of a many-body quantum system whose 
dynamics is described by a master equation 



d t p = Sp. 



(2) 



This equation describes the dynamics of an open system 
coupled to a Markovian reservoir, where p is the system's 
density operator. The superoperator <S contains two 
parts: one is related to the system Hamiltonian (even- 
tually rcnormalizcd due to the interaction with the envi- 



ronment) and the other with the dissipation induced by 
the environment. Under the appropriate conditions, the 
system evolves to a steady state p ss , which corresponds 
to a (right) eigenstate of <S with eigenvalue 0. Note that 
this eigenvalue may be degenerate, or there may be other 
eigenvalues with zero real part, in which case the steady 
state is not unique. In case this does not happen, the 
steady state is unique. Then, the other eigenvalues A 
of <S have a negative real part, and the smallest abso- 
lute value of them, A, determines the asymptotic decay 
rate (ADR), that is, the rate at which the steady state 
is reached. A phase transition in the steady state, where 
its properties abruptly change when one slightly changes 
a parameter in the master equation will be accompanied 
by the vanishing of A. This is typically referred to as a 
"dissipative quantum phase transition" (see, for example, 
3, 5, 17, 19, 20]). There is a natural analogy between dis- 
sipative and (closed-system) quantum phase transitions: 
A unique ground state of the Hamiltonian is analogous 
to a unique steady state. The appearance of a phase 
transition is signaled by the vanishing of the gap or A, 
respectively. 

Apart from its role in reflecting the appearance of a 
phase transition, the quantity A can play an additional 
role. It also represents a physical property of the sys- 
tem, namely the rate at which the steady state is ap- 
proached asymptotically or its response to perturbations 
in the steady state. This quantity may change abruptly 
itself. In that case, we can talk about a dynamical transi- 
tion, since a small change in the system parameters may 
lead to an abrupt change of the dynamics of the system. 
Actually, such a transition may in principle occur even if 
A remains finite, and thus it is a different property than 
the transitions generally studied in this context. 

In this paper we investigate both kinds of transitions 
for simple fermionic systems. We concentrate on mas- 
ter equations in which the Hamiltonian part is at most 
quadratic in fermionic creation and annihilation opera- 
tors. Additionally, we consider two kind of dissipative 
parts in terms of their dependence on such operators: (i) 
general quadratic and (ii) quartic, but with some condi- 
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tions (in particular, that they correspond to Hcrmitian 
Lindblad operators). In the first case, the dynamics can 
be exactly solved |2ll - l24l ] which has been exploited in 
several recent works to study the interplay of dissipa- 
tion and critical Hamiltonians in Id fermionic systems 
[ill [H, [25[. In the second case, even though the full 
dynamics cannot be obtained, we will show that it is 
nevertheless possible to exactly determine the dynam- 
ics of certain expectation values, from which dynamical 
and steady state properties can be obtained. In this last 
case we will present analytical examples where dynamical 
transitions occur [26| . This situation has also been stud- 
ied in j27l - [29j with particular regard to transport through 
a dephasing spin chain, where exact solutions of the as- 
sociated master equation could be obtained. 

Even though the formalism we develop is relatively 
general, we will illustrate it with explict examples. In 
particular, we consider Hamiltonians which are inti- 
mately connected to physical situations that can be ob- 
tained in the lab, namely anisotropic XY spin chains 
in transverse magnetic fields, and that a Jordan- Wigner 
transformation maps it into a fermionic Hamiltonian. 
This family of Hamiltonians displays the prototype of 
a continuous phase transition [l8[ . The dissipative terms 
we consider can also be understood as particular physical 
processes occuring in the spin chain through its interac- 
tion with an environment. Note that our framework also 
applies to the systems studied in [U [22|, [2?], I2II and 
for the quadratic dissipative terms is related to |l9l. l23j , 
where generic dissipative phase transitions are analyzed. 

This paper is structured as follows. In Sec. |H] we in- 
troduce the Lindblad master equation which allows to 
describe decoherence due to the weak interaction with a 
Markovian bath and present the covariance matrix for- 
malism, which allows the exact treatment of quadratic 
fermionic systems. In Scc. lIIII we extend this formalism to 
decoherent systems with linear and Hermitian quadratic 
Lindblad operators. Then we come to the calculation 
of the steady states and the ADRs for relevant interest- 
ing examples in this framework in Sees. IIV1 fVl and IVII 
Here we explicitly demonstrate the presence of dissipa- 
tive phase transitions. In Sec. IVIII we propose a possi- 
ble implementation with cold ions before concluding in 
Sec.EnU 



II. NOTATION AND METHODS 



In this section we introduce our tools and notation, 
namely the Lindblad master equation and the fermionic 
covariance matrix formalism which is ideally suited for 
describing quasi- free fermionic systems (see Sec. Ill Cl) . 



A. Lindblad Master Equation 

We treat the decoherence of a system in contact with 
a heat bath with the Lindblad master equation (30j 

d t p = Sp 

= -\ [ H "°] + E ( L > LQt - \ { LQtLa ^}) . ( 3 ) 

a \ / 

where p is the density matrix of the system, H is its 
Hamiltonian, and the Lindblad operators L Q determine 
the interaction between the system and the bath. This 
dynamical equation for an open system can be derived 
from two different points of view [3l| : First, it can be 
derived from the full dynamics of system and bath. Here 
three major approximation have to be used: The states 
of system and environment are initially uncorrelated, 
the coupling between system and bath is weak (Born 
approximation), and the environment equilibrates fast 
(Markov approximation). Second, any time-evolution 
given by a quantum dynamical semigroup (i.e., a family 
of completely positive, trace preserving maps e t , which is 
strongly continuous and satisfies ete s = et+ s ) is generated 
by a equation of the form Eq. [3J 

We characterize the decoherence dynamics with the 
steady state and the asymptotic decay rate. A steady 
state density matrix po of the master equation © fulfills 

d iPo =Sp = (4) 

and is the (generically unique) eigenvector with eigen- 
value of the Liouvillian superoperator <S. The approach 
to the steady state is then governed by the non-zero 
eigenvalues (and eigenvectors) of S, all of which have 
non-positive real part for Liouvillians of Lindblad form. 
Of particular interest is the eigenvalue with the largest 
real part, since it governs the long-term dynamics. We 
refer to the absolute value of this largest real part as the 
asymptotic decay rate (ADR) and denote it by A: 

A(«S) = max{|ReA| ± : 3p x : 5(p A ) = Xp x }. (5) 



B. Quasifree Fermions and Spins 

We consider systems with N fermionic modes de- 
scribed by creation and annihilation operators and Oj . 
These operators obey the canonical anti-commutation re- 
lations 

{aj,a k } = 0, {a],a k } = 5 jk . (6) 

Equivalcntly, we can use Hcrmitian fermionic Majorana 
operators 

c j ,o = a] + a jl c jA = (-i)(a]-a j ), (7) 
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which as generators of the Clifford algebra satisfy the 
anti-commutation relations 

{cj,u> c k,v} = 2SjkS uv . (8) 

We consider fermionic Hamiltonians that are quadratic 
in the Majorana operators. They describe quasifree 
fermions and are known to be exactly solvable. We pa- 
rameterize them with the real antisymmetric matrix H 

H ^ t Hjk,uvCj,uCk.v- (9) 

jkuv 

The 2x2 matrix Hjk = {Hjk tUV ) uv describes the coupling 
between the modes j and k. 

All eigenstates and thermal states of such a quadratic 
fermionic Hamiltonian are Gaussian, i.e., they have a 
density operator which is the exponential of a quadratic 
form in the Majorana operators. Gaussian states remain 
Gaussian under the evolution with quadratic Hamiltoni- 
ans. 

In the following, we will mostly concerned with trans- 
lationally invariant systems and nearest-neighbour inter- 
actions. In terms of the matrix H the former means that 
Hjk depends only on the difference j — k and we write 
for short 

H ]k = H^ k , (10) 

while the latter implies that H s = for s > 1. We work 
with periodic boundary conditions, so j — k is understood 
modulo TV. 

An important reason to study one-dimensional 
fermionic systems with quadratic Hamiltonian is their 
intimate relation to certain types of spin chains: The 
Jordan- Wigncr transformation [321 ] maps fermionic oper- 
ators onto Pauli spin operators via 

j'-i i-i 

k=l k=l 

Under this transformation some spin chains are mapped 
to spinless quasifree fermionic systems which can be 
solved exactly. A prominent example is the anisotropic 
XY chain in a transverse magnetic field [l8| with the 
Hamiltonian 

N 

H = - J E [C 1 + y)44 +1 + (! - y>>i +1 ] 

3 = 1 

N 

+ (12) 

3=1 

where B is the magnetic field, J the ferromagnetic cou- 
pling, and 7 the anisotropy parameter. Closed systems 
governed by this Hamiltonian show a quantum phase 
transition at B = 2J in the thermodynamic limit and 



the behavior in the presence of dissipation is studied in 
Sec. lVTBl 

We are interested in dissipative (open) fermionic sys- 
tems, with dynamics described by a Lindblad master 
equation, characterized by a set of Lindblad operators 
L a . We consider two classes of Lindblad operators: 
firstly, those given by arbitrary linear combinations of 
the Majorana operators (linear Lindblad operators) 

= / ] L™ u Cj )U , i"„ £ C, (13) 

ju 

and secondly, those represented by quadratic expressions 
in the Majorana operators which are in addition Hcrmi- 
tian (Hermitian quadratic Lindblad operators) 

L " = I E L %,uvC3,uCk,v (14) 
jkuv 

with the real and antisymmetric matrix L a . 



C. Covariance Matrix Formalism 

Now we present a framework in which the dissipative 
dynamics of the Lindblad master equation (|3|) can be 
solved exactly. 

For every state of a fermionic system, its real and an- 
tisymmetric covariance matrix (CM) is defined by 

Tjk,uv = tr (^p- [cj tU , Ck, v ]^J ■ (15) 

The magnitudes of the imaginary eigenvalues of F are 
smaller than or equal to unity (r 2 < — 1). 

For Gaussian states the correlation functions of all or- 
ders are related to the covariance matrix through Wick's 
theorem (33j . In particular, pure Gaussian states p = 
\^>)( 1 $\ satisfy r 2 = — 1. In our notation Tj k denotes a 
2x2 matrix that describes the covariances between sites 
j and k. 



III. LINDBLAD MASTER EQUATION IN THE 
COVARIANCE MATRIX FORMALISM 

The covariance matrix formalism is especially useful 
if the operative dynamics leads to closed equations for 
the CM, which is the case for the two kinds of Lindblad 
operators Eqs. (|13I14[) that we study in the following. 



A. Linear Lindblad operators 

We consider a system with quadratic Hamiltonian 
given by the antisymmetric matrix H [cf. Eq. (|9])] and 
linear Lindblad operators as defined in Eq. ()13l) . Using 
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the anti-commutation relations (J5J we determine the dy- 
namical equation for the covariance matrix L from Eq. ([3]) 
and obtain: 

d t T= [H,T] -^]{|L a )(L"| + \L a *)(L a *\ 7 T} 

a 

-2i(\L a )(L a \ - \L a *)(L a *\) , (16) 

where \L a ) denotes the vector formed by the coefficients 
L" u in Eq. (fT3|) and \L a *) its complex conjugate. In 
terms of |r), the vector of components of T, this equation 
becomes 

d t \T) = S\T) - \V) = (H-M) |r> - |V), (17) 
with the superoperators 

H= (H®1-1®H T ), (18) 

m = Yj (l L ")( LQ l ® 1 + 1 ® (l iQ )( iQ D T + c - c -) . 

a 

(19) 

|V) = 2^(|L")®|L Q )-c.c). (20) 

a 

Note that % is anti-Hermitian and M is Hermitian and 
positive semi-definite. The steady state covariance ma- 
trix [see Eq. fl3)] satisfies 

(H-M)\T ) = \V). (21) 

Deviations 1 517) = |r) — |r ) then obey 

d t \6T) = {U-M)\5T) (22) 

and the approach to the steady state is governed by the 
the right eigenvalues of the superoperator S = % — A4, 
satisfying 

S\Ti) =Xi\Ti). (23) 

The eigenvalues whose real parts are closest to zero thus 
determine the asymptotics of the decoherence process. In 
the following, we refer to 

A = max{|ReAi| ^ : 3L, s.th. (S - X l )\T l ) = 0} , 

(24) 

i.e., the asymptotic decay rate on the level of covariance 
matrices simply as ADR. 

B. Quadratic and Hermitian Lindblad operators 

The second class of master equations leading to closed 
equations for the covariance matrix is of the form Eq. §3§ 
with Lindblad operators that are quadratic and Hermi- 
tian, as in Eq. (|14j) . Lindblad equations with Hermitian 
Lindblad operators describe the dynamics of systems in 
contact with a classical bath. Let us choose a fluctuating 



external field as the source of decoherence (see Sec. IVIII) . 
If, additionally, the Lindblad operators are quadratic, the 
fluctuating Hamiltonian is quadratic. Thus in this case 
Gaussian states evolve into mixtures of Gaussian states 
under such evolutions and we can expect a closed equa- 
tion for the covariance matrix. 

Before discussing the master equation in the CM for- 
malism, let us first determine in general the steady state 
density matrices [see Eq. of a master equation with 
only Hermitian Lindblad operators. In that case, we can 
rewrite the master equation in terms of \p), the vector of 
components of p as 

d t \p)=s\p)= ^-^E( £Q ) 2 ) \p)- ( 25 ) 

with the superoperators 

H = -i (H ® 1 - 1 $ H T ) , (26) 
C a = L Q <S>l-l<S>L aT . (27) 

We observe that the superoperator ti is anti-Hermitian 
and that the superoperators C a arc Hermitian, so that 
the (C a ) arc Hermitian and non-negative. 

We consider all complex valued vectors |p) instead of 
just the ones corresponding to positive density matrices 
with trace one. Therefore, we have to check after the 
calculation if our results correspond to physically mean- 
ingful states. The steady states satisfy 

(/'o|(W-^E( £Q ) 2 )l' o)= - ( 28 ) 

a 

As stated above, H is anti-Hermitian and all (jC q ) 2 are 
Hermitian. Applying these properties we can conclude 
from Eq. ([25]) that 

( Po |E(£ a )>o) Hpo|«|a)> =0 (29) 

ct 

holds. It follows from the non-negativity of (£ a ) 2 that 

(r)>„) = ov a . (30) 

Because the C a can be diagonalized this implies C a \po) = 
0. It follows that n\po) vanishes identically. In terms 
of matrices po ; we can summarize these conditions for 
steady states 

[H, Po ]-[L Q ,po] = Va. (31) 

It can be verified with Eq. ([3]) that this condition for 
steady states is not only necessary but also sufficient. To 
summarize, steady states for Hermitian Lindblad opera- 
tors correspond to density matrices commuting with the 
Hamiltonian and all Lindblad operators. Therefore, they 
are the identity up to symmetries shared by the Hamil- 
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tonian and the Lindblad operators. 

Let us now return to exactly solvable systems in the 
CM formalism. For quadratic and Hcrmitian Lindblad 
operators and quadratic Hamiltonians the Master Equa- 
tion ([3]) becomes 



(32) 



We can again reformulate this equation for the vector of 
components T) 



9 t |r> = s|r> = ^-^E( £Q ) 2 ) 



r), (33) 



with H as in Eq. ([THJ and C a = L a ® 1 - 1 ® L a . 

Since we found that steady states are trivial for Hcrmi- 
tian Lindblad operators, we concentrate on the asymp- 
totics of the decoherence process. It is studied through 
the eigenvalues A; of the superoperator S, and in partic- 
ular its ADR as defined in Eq. (j24|) . 



C. Translationally invariant Hamiltonians 

Naturally, translationally invariant systems are best 
treated in a Fourier transformed picture. Any real an- 
tisymmetric matrix can be transformed into a real and 
antisymmetric block-diagonal matrix by an orthogonal 
transformation O. For the Hamiltonian matrix H this 
means 



TTl 

mn.uv 



(OHO 7 



e rn 



(34) 

where the real number e m are the energies of the elemen- 
tary excitations. We, however, transform the Hamilto- 
nian matrix with the unitary Fourier transform 



[UHU 



1 2xi 

V N 



(35) 

The resulting matrix H is anti-Hcrmitian, but not real. 
For translationally invariant systems, for which the 2x2 
matrices Hjk in Eq. §§§ depend only on j — k, the matrix 
H is block-diagonal with 



N-l 



e N 



(36) 



s=0 



The block-diagonal is parameterized according to 

H nn = y , k n ,l n e R, h n e C. (37) 

For later use, we observe the properties 



u - h* k — —h 1 - —I 



(38) 



which follow directly from Eq. (f3"6")l for real H s . 

For a system that is also invariant under reflections (in 
real space) H s = — if J holds (in addition to ff_ s = — Hf 
implied by antisymmetry). In that case, we have H nn = 



— u n 



In = 0. 



(39) 



The spectrum of the Hamiltonian matrix determines 
the elementary excitation energies 



(40) 



It will be necessary to transform the covariance matrix T 
accordingly, defining 



r = uruK 

By minimizing the energy expectation value 
(E) = Tr (H T T) = Trf/f+f), 



(41) 



(42) 



we find the covariance matrix for the ground state. In 
the case k n l n < \h n \ 2 it is 



r u - s 

1 mn — u mn 



-V2 % 



and otherwise 

r"„„ = -i<5 ln „sign (k„ + l n ) 1 2 . 



(43) 



(44) 



For translationally invariant and reflection symmetric 
systems k n l n = holds, thus k n l n < \h n \ 2 is fulfilled 
in such systems. Since the XY chain Eq. [T^is reflection 
symmetric, we can concentrate on the case of Eq. (|43|) . 
Specifically, we obtain for the Hamiltonian Eq. [12] that 



h„ = -2B + 2 J 



(l + 7)eTT" + (l- 7 ) e - 

h — i — n 



. (45) 
(46) 



which contains a continuous quantum phase transition at 
B = 2 J, where the gap closes and an elementary excita- 
tion energy e„ = \h„\ = exists. This Hamiltonian will 
be further discussed in Sec. [VI] 



IV. LINEAR LINDBLAD OPERATORS 

Now we apply the formalism introduced in the previ- 
ous Sections to some simple cases of physical interest. 
Here we choose the simplest examples, i.e., linear Lind- 
blad operators (see Sec. MI Af . We study two settings. In 
Sec. IIV Al we look at systems without any unitary evo- 
lution, observing dynamic transitions when tuning the 
strength of competing decoherence processes. Here we 
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enrich our presentation with an example for dissipativc 
state engineering. In Scc. lIVBl wc consider open systems 
governed by a Hamiltonian, which describes a quantum 
phase transition itself, and show that the dissipative sys- 
tem undergoes a transition for the same values of the 
system parameters. 



the Fourier transform ()35l) 

d t Y = -g 2 (p 2 +v 2 )Y 

( N 



9 2 V 



@cos(27m/iV)cr z ,r 



N 



?7 = 1 

N 

2g 2 Hi'(^ism(2irn/N)<T :l 



(51) 



A. Purely dissipative systems 



In this case, a spectrum of decoherence rates g 2 {^ 2 + v 2 ± 
2//i/[cos + cos 2*p] i ^2 + y i ± 2 ^[cos 2=» _ cos 2^p j} 

exists with a "gap" g 2 (fi — is) 2 . The unique steady state 
is 



The simplest example of two competing decoherence 
processes generated by linear Lindblad operators is 



gfj,a 



LI 



(47) 



acting on site a € {1, ... ,N}. It describes the compe- 
tition between particle-loss and particle-gain processes. 
We observe that the Master Equation (fT5| without the 
Hamiltonian (H = 0) is diagonal in real space 



N 



d t r = -.9 V + ^ 2 )r - 5 V - v 2 ) ©(ia,,). (48) 



In this simple case the master equation is already diag- 
onal and we read off the single decoherence rate A = 
g 2 (/J, 2 + v 2 )- Solving the master equation for c^To = 
gives the unique steady state covariancc matrix 



r n = - 



A t 2 - v 2 

fi 2 + IS 2 



N 





1 

-1 



(49) 



which is block diagonal. This state is characterized by 
the particle number {a* a a a ) = v 2 /(^ 2 + is 2 ) at all sites. 
For pure particle- loss processes (v = 0), all sites are un- 
occupied (a^a a ) = in the steady state, while for pure 
particle-gain processes (p = 0), all sites are occupied 
(a£,a Q ) = 1. At fj, = v the steady state is the unpolar- 
ized completely mixed state. Not surprisingly, the system 
does not display any phase transition. 



r 



/i 2 + z/ 2 



N 

0^ 



2/i v 



N 



V 9 
[I 1 



0isin(27rn/iV)a- :j 



(52) 



This state is a paired fermionic state according to the 
definition of Kraus et al. (34|- Paired states show two- 
particle quantum correlations that can not be be repro- 
duced by separable states (mixtures of Slater determi- 
nants). It is proven in [34[ that Gaussian states are 
paired iff Qki = (%[a,k,ai]) 7^ 0. This condition expresses 
the fact that separable states are convex combinations of 
states with a fixed particle number. For the covariance 
matrix (|52l) we get 



!ki = 



1 /_iz> • sigii(/V 

2 fj,*+v* 





if \k-l\ = 1 

if \k-l\ ^ 1 



(53) 



We conclude that ([50)) generates paired states, except 
for the trivial cases [i — or v — 0. Note that even 
though the gap closes at fi = v (where maximal pairing 
is created) there is no phase transition at this point. 



B. Dissipative systems with Hamiltonians 



More interesting may be the case in which dissipation 
can also induce correlations. A simple example of this 
kind is provided by the Lindblad operators 

L" = g ( M a a + (50) 

acting on nearest neighbors. This set of Lindblad opera- 
tors generates a master equation, which is diagonal after 



A different form of transitions can arise in the pres- 
ence of a Hamiltonian when tuning the parameters of the 
Hamiltonian. To show this, we solve the evolution of the 
Lindblad master equation (|16p with a general quadratic 
and translationally invariant Hamiltonian [see Eqs. Q 
and ((37]) ]. We choose the local Lindblad operators (|47|) . 
again because they are the simplest example. The diag- 
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onal master equation in Fourier space becomes 

N 



01 ik n h n 
\-h* iL 



r 



N 



g ^ + ^)T-g\^-u 2 )^ia y . (54) 



n=l 



The corresponding steady state covariance matrix in the 
weak-coupling limit g — > is [35| 



r 



JY 

e 



Re (h n ) 



H 2 + v 2 ^ 1 (k n -l n ) 2 /4+\h n \< 



i(k n In) 1 2 h n 

-h* n -i(k n - l n )/2 



(55) 



Transforming back to To [and using Eq. (|38[) ] we can read 
off the particle number (2ajaj — f ) = (To)jj,oi as 



(ro)jj,oi - o 



1 1? - v 2 1 



JV 

E 



Re (h n ) 



^ + U 2 N ^ (k n -i n y/A+\h n [< 



(56) 



Based on this result we can now discuss how non-analytic 
behavior in the steady state correlates with critical points 
of the system. A vanishing denominator in Eq. (|56p is not 
a priori a sufficient condition for non-analytic behavior 
because the numerator might vanish at the same point. 
This is relevant for interesting examples with k n — l n = 0, 
e.g., the XY chain in Eq. ([4"5"j) . We give a rigorous dis- 
cussion in the following. In the thermodynamic limit, 
the sums over expectation values in Eq. (|56[) can be re- 
placed by a loop integral around the origin of the complex 
plane with radius one, where the integration variable is 
z = cxp ( 2 j^nJ . This is possible because h n , k„, and l n 
are Fourier series. For local interactions, the denomina- 
tor of the integrand is a polynomial in z [see Eq. ([55]) ] 
and thus has a finite number of distinct roots. Apply- 
ing the residue theorem, a non-analyticity in (fljaj) is 
possible only if a residue of the integrand, i.e., a root of 
its denominator, moves through the integral contour in 
the complex plane as a function of some external pa- 
rameters. This happens for a vanishing denominator 
\h n \ 2 + (&n - ^i) 2 /4 = for some real n <E [0, N). In the 
special case of a reflection symmetric system fc„ + l n = 
this coincides with a vanishing energy gap e n = [see 
Eq. (|40p ]. a signature for a quantum phase transition. 
To summarize, for a reflection symmetric system with 
\h n \ 2 + (k n — l n ) 2 /4 = in the weak-coupling limit a 
quantum phase transition occurs in the dissipative sys- 
tem for the same parameter values as in the correspond- 
ing closed system and is signaled by a non-analyticity in 
(ajCij). This calculation is explicitly performed in section 
EE] for the XY chain 



V. QUADRATIC AND HERMITIAN LINDBLAD 
OPERATORS 



In this Section we turn to the dynamical properties of 
the Lindblad master equation with quadratic and Hcr- 
mitian Lindblad operators as introduced in Sec. IIIIBI 
In the study of closed systems, quantum phase transi- 
tions are signalled by non-analyticitics in ground state 
expectation values. In the dissipative case the steady 
state is the analog of the ground state. However, we 
have shown in Sec. MI Bl that in the case of Hermitian 
Lindblad operators the steady states are trivial and thus 
cannot evidence a phase transition. Therefore we turn 
to the ADR, which determines the long-time dynamics 
of the decoherence process. We identify non-analytical 
behavior of this rate both in the absence of any Hamilto- 
nian (see Sec. IV Ap for competing decoherence processes 
and for non-zero Hamiltonian, in which case phase tran- 
sitions of the corresponding closed system are reflected 
in a "dynamical transition" of this rate (see Sec. IV Bj) . 



A. Purely dissipative systems 

A particular simple set of local and quadratic Lindblad 
operators is 

2 

L z = [Ca,U c afi] , (57) 

1 

L £ = [ C «+!.0: C a,l] ■ (58) 

In this case the Lindblad equation ([52"]) becomes 

dtTkLuv = - 4g 2 ii 2 r kl , uv (1 - 5 U ) (59) 

— Ag 2 V 2 Ykl.uv(^ — ^2k+u+1.2l+vSk+l.l 

— <52fc+K-l,2Z+u<5fc-l,z), 

We can read off the decoherence rates — 4g 2 (/i 2 + v 2 ), 
-4g 2 n 2 , and -Ag 2 v 2 . Thus, the ADR 



A 



4g 2 /j, 2 if fi < v 
Ag 2 v 2 if v < (j, 



(60) 



undergoes a dynamical transition as a function of fi/v at 
fx = v. 



B. Dissipative systems with Hamiltonian 

Now we add a quadratic Hamiltonian and calculate the 
asymptotic decay rate A in the limit of small couplings 
to the environment g — > 0. First, we derive it for the 
quadratic Lindblad operators from Eqs. (|57l58p for v = 
and fi = 1 . Later we will present the results for the case of 
arbitrary /z and v. For translationally invariant systems 
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the Fourier transformed master equation (|59|) is 
d t f M = (ST) M = [H,f] kl 

- 4g 2 (f kl - 1 ^ f r A- s , fe -*J , (61) 

V r,s=l / 



with the unitarily transformated superoperator S to 
S= (U®U)S(U®U)\ 



(62) 



with U from Eq. ([55]) . For weak couplings between sys- 
tem and bath g — > 0, the eigenvalues of S (and thus of 
iS) can be determined by first order perturbation expan- 
sion. To this end we first diagonalize the unperturbed 
Hamiltonian part of S 



kl 



H kk r 



fcfcJ- kl 



XT 



kl- 



(63) 



where we use the notation introduced in Eq. ([3"rJ| for 
the Hamiltonian H . The 4N 2 eigenvalues \ mna (m, n = 



1,...,N, a = L. 


. , 4) are 






ymnl 






(64) 


yrnn2 


— ^ (j%m 0m "i" 




(65) 


yrnnZ 


= i (a m - a n + (3 m -f 


AO, 


(66) 




— ? (Orn — ftm ~ 


AO, 


(67) 


with 














(68) 


Pm = 


= \/|/lm| 2 + (fern ~ 


2/4. 


(69) 



The corresponding eigenmatrices are denoted as A mno 
with nonzero elements A£]™ a only for m = k and n = I, 
i.e., A™ a = S mk S n iA™™ a . Perturbation theory demands 
to calculate the matrix elements of the perturbative part 
of S, ~4g 2 (S mk S nl S ab ~ Ptfr/N) [see Eq. JH)], with 



Tjrn.na 
^klb 



-T. 

2 ^ 



N 



q,r,s,t— 1 



(£«)'|A W6 >+J\rW»J««6 

(A£T 



(70) 



Thus the eigenvalues of S are determined by those of the 
Hermitian matrix P and the largest eigenvalue smaller 
than N of P (restricted to a space of degenerate eigen- 
values X mna of [H, •]) determines the ADR. We denote 
it by Ap and thus have that the asymptotic rate is 



A = Vl - 



N 



To find Ap, note that the matrix 



elements of P fulfill \P^ a \ < 1. Thus an JV-fold degen- 
eracy of X mna is required for Ap = O(iV). Gencrically, 



this is possible only for the eigenvalue X mna = 0, i.e., 
m = n and a = 1,2. The corresponding eigenmatrices 
are 



\rnml _ r 



1 



mk"ml 



A iM — 0,nkO,nl 



J_ A 

V2 10 1 



2 

h* 



h m 



(71) 
(72) 



As the eigenmatrices A mm2 give eigenvalues equal to A" 
and only, have no overlap with physical covariance ma- 
trices, and yield P k l^ 2 = 0, we focus on the matrices 
A mml . The corresponding part of the perturbation ma- 
trix is 



p prom 

r mn — r nnl 



1h m h* n + 1h* m h n — (k m — l m j (k n — l n ) 

4/3 m /3„ ' 
(73) 

We diagonalize this matrix by introducing the three vec- 
tors la), |o), |c) e with the components 



Im( /i T - 



Re(/i m ) 



(74) 



^Pm Pm Pm 

and writing P mn in terms of these unnormalized vectors 

P = \c)(c\ + \b)(b\ - \a)(a\. (75) 

We now exploit the symmetries of h n , k n , and l n stated 
in Eq. (f3"5|) . First we observe that | c) is orthogonal to | a) 
and | b). We have chosen the covariance matrices corre- 
sponding to the three vectors (|74[) anti-Hermitian, since 
this matrix remains anti-Hermitian even in the complex 
vector space. After transforming back into real space the 
ones corresponding to | a) and | b) are purely imaginary so 
that they have no overlap with any physically meaningful 
real and antisymmetric covariance matrix. Only the ma- 
trix corresponding to c) is real and antisymmetric and 
given by 



JV 



y-v Rc[h ri 



Therefore, it determines the ADR. We get 



U f A nnl U. (76) 



A, 



JV 

E 



IM 2 + (fc m -U 2 /4' 



and thus 



4g 2 4Im(/i m ) + (fc m - / TO ) S 

rn—1 



as the general from of the ADR. 



/ ¥ 



(77) 
(78) 



We can extend our analysis to systems with the general 
Lindblad operators Eqs. (|57I58I) and find in an analog 
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way the two lowest decay rates 



ELi h* m lh m in the thermodynamic limit by in- 
troducing the complex variable z = exp ( — 1mm /N) 



,,2 



± 



(79) 



with 



El. V 



Re(hrr,) 2 



N ^ B 2 

m—l " 



v 2 ^ Re(/i m exp(— 2Trim/N)Y 



m—l 
N 



[iv ^--^ Re(ft, m )Re(ft, m exp(— 2nim/N)) 



N 



rn— 1 



/? 2 



(80) 



We can now argue that the ADR itself reflects the criti- 
cality of the system. The argument is completely analo- 
gous to the one given in Sec. IIVBI If the denominator be- 
comes zero, we can expect a non-analyticity expressions 
Eqs. (|80|) . In particular, in the reflection symmetric case 
k n + In = 0, where the denominator agrees with the ele- 
mentary excitation energies [e 2 = \h n \ 2 + (k n — /„) 2 /4 = 0, 
see Eq. (|40)) ] the non-analyticity in the ADR signals the 
presence of a quantum phase transition in the Hamilto- 
nian itself. 



VI. EXAMPLE HAMILTONIANS 

In this section we will revisit the results obtained for 
the steady state and the ADR for linear and quadratic 
Lindblad operators in Sees. II V Bl and fVBl for the specific 
Hamiltonian (|12p of the quantum XY chain. 

The energies of the elementary excitations of this 
Hamiltonian are e„ = \h n \. Thus, for the XY chain in 
Eq. (|4"5j) the gap closes at B = 2 J in the thermodynamic 
limit and the quantum XY chains exhibit a phase tran- 
sition at this point. In fact, these models constitute the 
archet ypa l example of a continuous quantum phase tran- 
sition [18l | . In this chapter we want to find properties of 
the dissipative dynamics signaling this phase transition. 



N^oo N h r 
m—l 



2ni 



dz 2J(1 - j)z 2 - 2Bz + 2.7(1 + 7) 
| =1 ~ 2J(1 + I )z 2 -2Bz + 2j(l -7) 



(82) 



where the integration contour is a circle of radius \z\ = 
1 around z = in the complex plane. The complex 
integrand is analytic except for three distinct poles at 



0, 



1 



B± v /S 2 -4J 2 (l-7 2 ) 



(83) 



2J(l+ 7 ) 

The contour integral is determined by the sum over the 




FIG. 1. The poles z°,z ± [see Eq. ^] are plotted for J = 
1,7 = ±0.1. As B is changed from to 20 the poles z + (z~) 
for positive anisotropy 7 = +0.1 move along the blue (red) 
solid curves and z + crosses the contour at the critical value 
B = 2J. For negative 7 = —0.1, z~ crosses at B = 2.1 . At 
the crossing the integral Eq. ()82|) changes non-analytically. 



A. Linear Lindblad operators 



Let us now apply the findings from Sec. IIVBI and Eq. 
([55]) to the example system defined in Eq. (|4"5j) which 
contains a quantum phase transition at B = 2J. Then 
the particle numbers become for k n = l n = 



(24 a n - 1) = - 




For 7 = we easily obtain A = (since by Eq. (|4"5)l 
h n is real in that case and then by Eq. ((78)) A is zero 
for k n = l n = ). For 7 7^ we evaluate the sum 



residues at those poles which are inside the contour (\z\ < 
1). z° is always inside this contour. In the case 7 > 0, 
z + is inside the contour for < B < 2J and outside for 
B > 2J, while z~ is always inside the contour. In the 
case 7 < 0, z~ is inside the contour for B > 2 J and 
outside for < B < 2J, while z + is always outside the 
contour. So residues cross the contour at the quantum 
phase transition B = 2 J (because then h n = for some 
n), leading to a non-analytical behavior in the particle 
density of the steady state. 

After applying the residue theorem we get the particle 
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FIG. 2. Asymptotic decay rate A [see Eq. J75J] of the XY 
chain (|12|l for different anisotropy parameters 7 as a function 
of the magnetic field in the limits N — > 00 and g — s- 0. A 
phase transition in A is visible at B — 2 J for 7 / 0. 



number of the steady state 



(2at a n - 1) = 



1 

1-7 2 



1 - 



v / i-(¥) 2 d-7 2 ) 



B <2J 
B > 2J 

(84) 



for all 7, which does not depend on the sign of 7. For 
B < 2J the particle number in the steady state does 
not vary with the magnetic field, while its magnitude ap- 
proaches (/i 2 — v 2 )/(n 2 +is 2 ) for large magnetic fields like 

~ (J/B) . To summarize, the steady state undergoes 
a dissipative phase transition at B = 2 J signaling the 
phase transition in the system. 



B. Quadratic and Hermitian Lindblad operators 

As an example we study the anisotropic XY chain in 
a transverse magnetic field with the Hamiltonian given 
in Eq. (|12p . This translationally invariant Hamiltonian 
is Jordan- Wigner transformed to a quadratic fermionic 
Hamiltonian with Hamiltonian matrix H given by 



-2B 
2B 



2J(l-7j 
2j(l+ 7 ) 



Hi = 

H-! = 



2j(l+ 7 ) 
-2j(l- 7 ) 



(85) 
(86) 
(87) 



After Fourier transforming [see Eq. ([55]) ] this Hamilto- 
nian matrix assumes the form given in Eq. (|37[) with pa- 
rameters h n ,k n , l n given by (|45ll . 

We now apply the results from Sec. IVBI to the Hamil- 
tonian Eq. ([55]) and the Lindblad operators 



After a brief discussion of the steady states and a deriva- 
tion of the ADR in the thermodynamic (N — > 00) and 
in weak coupling (g — > 0) limits, we present numerical 
results of the system dynamics for finite N and g and 
compare them with our analytic predictions. 

First, we discuss the steady states of these systems (see 
Sec. MI Bp . From Eq. (pTTj) we have concluded that the 
steady state density matrix is the identity up to symme- 
tries shared by the Lindblad operators and the Hamil- 
tonian. A rigorous derivation of the steady states for 
this example could start from the ansatz that the steady 
state density matrix is diagonal in the Fock basis, follow- 
ing from [<t",/?] = 0. Then the commutator [H,p] =0 
must be exploited to get the steady state. 

As the Lindblad operators correspond to local particle 
number operators, the important compatible symmetries 
for the XY chains are the parity V = a\ . . . , discrimi- 
nating between an odd and an even number of particles, 
and the total particle number N = (l + ^2 °z) /2- For 
truly asymmetric XY chains 7 7^ 0.5 the parity is the 
highest symmetry compatible with the Lindblad opera- 
tors. In these cases the steady state density matrix is 
given by the identity in the two sectors of even and odd 
parity, the relative weight of these sectors is determined 
by the initial state. For the symmetric chain 7 = 0.5, the 
steady state density matrix is the identity only in the 
sectors with a constant total number of particles. Thus 
for 7 ^ 0.5 the steady state magnetization is (er^) = 
regardless of the initial state, whereas the magnetization 
of the initial state is conserved for 7 = 0.5. 

Second, we calculate the ADR ([75|) for the XY chains 
with Eq. ([4"5"j) analogous to the integration in Sec. IVI Al 

After applying the residue theorem we get the ADR 



M 
1+I7I 



i-(¥) (1-7 2 ) 



-1/2 



- 1 



B < 2 J 
B > 2 J 



(89) 

for all 7 in the case ^ = 1 and v = 0. It does not depend 
on the sign of 7 and is shown in Fig. [2] for several values 
of 7 G [0, 1]. For B < 2 J the ADR does not vary with 
the magnetic field, while its magnitude decays to zero for 
large mag netic fields like - (J/B) . The same behavior 
was found for the variance of the particle number in these 
models in a previous work (37| . To summarize, the ADR 
undergoes a dissipative phase transition at B — 23 sig- 
naling the phase transition in the system. 

The final result for the ADR ([89]) is valid in the lim- 
its N —¥ 00 and g — > 0. In this section we perform a 
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FIG. 3. Asymptotic decay rate A [see Eq. lf?8])] of the XY 
chain (|12[) for different coupling strengths g, 7 = 1, and N = 
100 as a function of the magnetic field B. For g < 0.1 (J/h) ' 5 
the results agree with the limit of weak coupling g — ► [see 
Eq. (jUJ)]. 
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FIG. 4. Asymptotic decay rate A [see Eq. ((89)) ] of the 
XY chain ()12|) for different system sizes N and 7 = 1, 
g = 0.01 (J/K) ' 5 as a function of the magnetic field B. For 
N > 50 the thermodynamic limit is reached except for small 
variations at the phase transition B — 2J . 



numerical diagonalization of the Lindblad master equa- 
tion superoperator S to compare the analytic result with 
the values for finite N and g. Furthermore, we extract 
the ADR from a simulation of the system dynamics and 
compare it with our prediction. 

In Fig. |3] we present the ADR for finite coupling 
strengths g. For g 2 < O.OU/fi the result of perturba- 
tion theory is in excellent agreement with the numerical 
diagonalization of the Lindblad master equation superop- 
erator. Deviations are strongest at small magnetic fields. 
We show the ADR A for different finite system sizes in 
Fig. [TJ Even in small systems with N = 10 spins the 




Time t / g 

FIG. 5. Evolution of the magnetization (oi) in time starting 
from the system ground state of the XY chain (|12|) . for dif- 
ferent magnetic fields B, g = 0.01 (J/h) "', and 7 = 1. The 
magnetization decreases exponentially in time. 
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FIG. 6. The ADRs A [see Eq. (89])] of the XY chain for 
7 = 1 for g = 0.01 (J/h) ' 5 and g — > (result of perturbation 
theory) as a function of the magnetic field B are compared 
with the late-time decoherence rates extracted from Fig. [5] 
The agreement between the ADR and the late-time decoher- 
ence rate shows the validity of our calculations for finite times. 



same qualitative behavior is found as in thermodynamic 
limit, i.e., the ADR signals the quantum phase transition 
in the system at B = 2 J. However, finite values of g and 

lead to a smearing out of the phase transition. 

We have defined the ADR through a diagonalization of 
the master equation, trying to describe the long-time dy- 
namics of the system. To demonstrate the deep relation 
between A and the dissipative dynamics, we extract the 
decoherence rate from a dynamical calculation (see Fig. 
[5|. Here we start from the ground state of the system 
and study the decay of the magnetization in time af- 
ter the system is brought into contact with a Markovian 
bath. In this example the exponential decay expected 
after long evolution times is nicely visible. In Fig. |5] we 
compare the extracted decay rates for different magnetic 
fields with the result of the diagonalization. We find an 
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exact agreement with the asymptotic decay rate numer- 
ically calculated with the same finite parameters. 

We can calculate the ADR for the XY chain for general 
values of [i and v in a similar way. In the spin picture 
the Lindblad operators are 



= 9^a 



(90) 

L " = 9 v(jX a °l+\ = [ Cq ' ' Ca + 1 - 1 ^ ■ ( 91 ) 
We find for the constants in Eq. (fT9"|) in the case 7 = 1 

tx/y 2 = 



B < 2 J 
B > 2 J, 



x 2 \2J, 



e/ (nu) 



B <2J 
B > 2J, 

B < 2 J 



(92) 
(93) 



D_ 

13 

^ B > 2J. 



(94) 
(95) 



In the symmetric case fi = is, the ADR A is constant 
A = — 4g 2 /i 2 . However, the next larger decoherence rate 
changes non-analytically: 



A_ = 



-2g 2 (i 2 ( 3 
-2 5 V (3 



if S < 2 J 
if B > 2 J. 



(96) 



pendent on their internal states. This state-dependent 
Coulomb repulsion can be designed to give the required 
spin-spin interaction. The spin state can be measured by 
fluorescence imaging of the ions. 

In this way the quantum Ising chain (39l . [loT ] and frus- 
trated Ising models [4l[ have been realized in recent ex- 
periments. In these experiments the ions were first cooled 
to their zero-point motional ground state and optically 
pumped into a certain spin configuration representing 
the ground state of the system without spin-spin interac- 
tions. Then the spin-spin interactions were adiabatically 
increased such that the system underwent a phase tran- 
sition. Finally, it was checked that the final state rep- 
resented the ground state of the simulated Hamiltonian. 
A large non-critical 2d Ising system has been simulated 
with ions in a Penning trap [421 ]. In the digital approach 
to quantum simulation with trapped ions, the elements of 
a general toolbox including Hamiltonian and dissipativc 
dynamics have been demonstrated [43|, EH . 

We describe in the following how to extend analog 
quantum simulation to include an incoherent evolution. 
The Lindblad master equation ([3]) with Hcrmitian Lind- 
blad operators L Q = ga^ (see Sec. IV B[) can be real- 
ized by introducing fluctuations of the simulated mag- 
netic field B a (t) = B a + SB a (t) [H as shown in the 
following. The local magnetic fields 5B a (t) should be 
uncorrelated between different sites 5B a (ti)5BP fo) = 
8 a p8B a (ti)8B a (t2). We restrict our derivation to a sin- 
gle Lindblad operator without loss of generality. Let, for 
example, 8B(t) constitute a Gaussian stochastic process 
of zero mean 5B(t) = with the time-correlations 



VII. EXPERIMENTAL REALIZATION 



6B(ti)6B(H 



SB 2 for \h-t 2 \<T 







for \h-t2\ > T. 



(97) 



We now discuss an experiment suited for the measure- 
ment of the asymptotic decoherence rate in spin sys- 
tems. The quantum simulation of spin systems with 
trapped ions was proposed in [38[ , where the spin degree 
of freedom is represented by two hyperfine levels. The 
magnetic field can be simulated either by directly driv- 
ing Rabi oscillations of the hyperfine transition or with 
position-independent Raman transitions induced by suit- 
ably aligned lasers. The spin-spin interaction is mediated 
via motional degrees of freedoms. State-dependent opti- 
cal dipolc forces (compare with state-dependent optical 
lattices) are generated by coupling the two hyperfine lev- 
els to electronically excited states with off-resonant laser 
beams. These dipole forces change the distance and con- 
sequently the Coulomb repulsion between two ions de- 



The correlation time T has to be much shorter than ev- 
ery process in the system (Markovian limit), i.e., H^HT < 
ujT <C 1, with the spectral width to of the Hamiltonian 
(difference between largest and smallest eigenvalue) and 
the superoperator % from Eq. (|2"6"|) . The averaged den- 
sity matrix evolves like \p(t)) =U (t) | /o(0)) , where the bar 
denotes the statistical average over the fluctu ating mag- 
netic field. The time evolution operator U (t) consists of 
contributions from % and 



Vi 



1 <g> V J 



(98) 



with V = a z . We can evaluate the statistical average of 
the time evolution operator in the interaction picture for 
the superoperators 
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U (f) = e w Texp (j dre-« 7 'V(T)e« 

OO 



n=0. 



tfti . . . dt„e-«'i V(ii)e"*i • • • e-« t «V(t„)e« t " 



E 

m— 



*>*i>— >*»>0 
l 2 "' I 



/ 



dt 1 ...dt rn e- Htl V' 1 e Htl 



U (i) = cxp m - - 



t>ti>— >t m >o 
l&B 2 ! 1 . 



2 ft 2 



-V 2 £ 



(99) 



with the time ordering operator 1~ . Between the sec- 
ond and the third line, we keep only even summation 
indices m = 2n (zero mean Gaussian process), eval- 
uate the statistical average at adjacent times tin—\ ~ 
tin < T (correlation time T), and neglect the terms 
exp (Hfan—i ~ t2n)) *C 1 (Markovian limit). In sum- 
mary, we have shown that the described fluctuations 
of the magnetic field generate Markovian dynamics [see 
Eq. (f2"5|) ] with Lindblad operators L Q = ga" = <?V and 
decoherence strength 
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SB 2 T 



(100) 



In the case of the anisotropic XY chain [see Eq. (|12p], the 
correlation time T is bounded by the width of the single 
particle excitation spectrum T _1 ^> max (4B/h, 8 J/h). 
In the recent experiment |39j 2 J/h ss B/h = 2tt x 4.4 kHz 
was used, but experimentally available laser intensities 
allow 2J/h ps B/h ps 2ir x 40 kHz. We propose to 
create fluctuations of the magnetic field with frequency 
T- 1 = 2ttx1.6MHz and variance JW/h 2 = {Q.2B/K) 2 
(2ir x 8 kHz) 2 . This would result in the decoherence 
strength g 2 pa 2 • 10 -3 J/fi and would require coherence 
times of order 2w/g 2 ps 25 ms. These coherence times 
can in principle be achieved in systems of trapped ions 

H. 



VIII. CONCLUSION 

We have investigated the dynamics of open quantum 
systems with regard to their steady states and asymp- 
totic decay. We have shown that insight into different 
phases can be gained by spectral analysis of the Liou- 
villian in analogy to how the spectrum of the Hamilto- 
nian reveals critical behavior in zero-temperature quan- 
tum phase transitions. 

To illustrate this point we have analyzed in detail the 



Liouvillian of open fcrmionic systems under a transla- 
tionally invariant, quadratic Hamiltonian, coupled to a 
Markovian bath. We treat master equations with linear 
or quadratic and Hermitian Lindblad operators. In both 
cases, the master equation leads to a closed equation for 
the covariance matrix from which the steady state co- 
variance matrix and the rates at which it is approached 
can be obtained exactly (see also [l9[ for an elegant and 
comprehensive treatment of both fermionic and bosonic 
linear open systems and their critical properties and [28[ 
for a detailed study of transport in spin chains under 
dissipation and dephasing). These results apply as well 
to a large class of Id spin systems that can be mapped 
to quasifree fermions by a Jordan- Wigncr transforma- 
tion. We have proposed an experimental realization of 
this quantum simulation with trapped ions. Numerical 
calculations show that our results for the weak decoher- 
ence limit do apply to such finite systems. 

We focused on the limit of weak decoherence (g — > 0) 
and showed how to deduce information about critical 
points from the spectrum of the Liouvillian. In particu- 
lar, the ADR A, i.e., the smallest non-zero eigenvalue of 
the Liouvillian, can serve an an indicator of phase transi- 
tions even if the steady state of the system is trivial and 
steady state expectation values thus cannot yield such 
information (as in the case of Hermitian Lindblad opera- 
tors). Depending on the decoherence process considered, 
the critical point can be reflected in the spectrum of the 
system's Liouvillian in the form of a closing gap (A —> 0), 
a degeneracy of A or non-analytic behavior of A. These 
results are summarized in Table HI 

With this work we suggest the possibility of detecting 
certain system properties through an observation of the 
decoherent dynamics: phase transitions in closed systems 
can be reflected in non-analytic changes of the ADR [25|, 
I2H |28| . More generally, since the ADR and other decay 
rates represent physical properties of the system, such 
non-analyticities can be seen as signature of a transition 
to a different dynamical regime. This suggests to study 
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TABLE I. Different dissipative systems studied, characterized 
by their Lindblad operators and Hamiltonian H . Relevant 
properties of the ADR A and the steady state are listed. x c 
denotes critical points of the Hamiltonian H . 



Lindblad Op 


ADR and Gap 


Steady State 


Hamiltonian H = 


fia a ,vdl 


gapped, no p.t. 


thermal 


fia a + va' a+1 


gap closes @ ^ = v, 
no p.t. 


paired 


i/i r i 


degenerate @ ^ = v 


oc 1 


^[Ca+1,0, c a,l] 






Hamiltonian H ^ 0: transl. invariant, critical at x c 




degenerate @ x c 


(era) non- 
analytic @x c 


i/t r l 


non-analytic @ x c 


oc 1 


^[Ca+1,0: c Q,l] 







the phase diagram of steady state correlation functions 
(A(t)B(t')) which will reflect these dynamical transitions. 
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